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Abstract 

We have developed an algorithm to imprint quantum fluctuations onto the initial transverse en- 
ergy density profile according to a given two-point covariance function. Using as an example 
MC-KLN initial conditions with added fluctuations satisfying the covariance function derived in 
[ 1 1, we find that effects from sub-nucleonic gluon field fluctuations on the eccentricity harmonics 
e„ vary strongly with the gluonic correlation length controlled by the saturation momentum Q s . 
Varying Q s over the range probed in Au+Au collisions at RHIC, we find gluon fluctuation in- 
duced enhancements of the eccentricity coefficients ranging from 10 to 20% in central collisions. 



1. Event-by-event fluctuations 

The importance of event-by-event fluctuations in the initial conditions context of heavy-ion 
collisions was first pointed out by Miller and Snellings who added nucleon position fluctuations 
to existing smooth (i.e. ensemble-averaged) Glauber initial conditions ||2][3). These fluctuations 
explain the experimental observation of non-vanishing anisotropic flow in central Cu-Cu and 
Au-Au collisions [4| and of odd flow harmonics {V3, vg, ...} |5|. 

Recently attention has turned to a new source of event-by-event fluctuations, namely fluctu- 
ations in the transverse distribution of color charge within the colliding nucleons. These fluctua- 
tions are evidenced by large multiplicity fluctuations observed in minimum bias pp collisions and 
suggest that sub-nucleonic degrees of freedom may also play an important role in determining 
the event-by-event geometry of the initial state in nucleus-nucleus collisions SHOD- 

2. From nucleonic to sub-nucleonic fluctuations 

In this work, we develop a toy model for imprinting sub-nucleonic fluctuations on the trans- 
verse energy density profiles produced in relativistic heavy-ion collisions. The study is motivated 
by recent work of Muller and Schafer in which they calculate the mean normalized covariance 
function Cov[e(r)/eo] for the transverse energy density fluctuations of gluon fields in central 
Au-Au collisions at 200 GeV Q. The authors approximate the collision system by two infinite 
slabs of nuclear matter with fixed gluon saturation momentum Q s for which they take the value 
corresponding to the nuclear thickness function in the center of a central Au-Au collision. 

We texture a given transverse energy density profile with additional fluctuations using a Turn- 
ing Band Gaussian random field simulator (TBSIM) |9 1 that includes several configurable covari- 
ance functions. The Miiller-Schafer covariance Cov[e(r)/eo] = (Ae(r)/eo) 2 is well described by 
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the Cauchy covariance included in TBSIM, 

Cov[r] = C(l +(r/a) 2 r 



(1) 



with fit paramaters C = 0.4679, a = 0.2878 fm and b = 1.7732. 

Using TBSIM we generate a large Gaussian random field (GRF) with the desired covariance 
on a 4000 x 4000 lattice with grid spacing Ax = 0.1 fm. The left panel in Fig. [T] shows a small 
slice of this GRF; in the right panel its two-point correlation function is compared to the Miiller- 
Schafer covariance and its Cauchy fit. 
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Figure 1: (Color online) Left: lOfmxlOfm section of the TBSIM Gaussian Random Field for the normalized energy 
density e(r)/eo~l with zero mean and a Cauchy covariance fitted to the Muller-Schafer covariance. Right: Two-point 
covariance of the GRF as a function of distance r at various azimuthal angles (symbols), compared with the Mtiller- 
Schafer function (solid) and its Cauchy Covariance fit (dotted). 



Due to the tail of the Gaussian distribution, there is a 7.2% chance that fluctuations about the 
mean fall into the unphysical region of negative energy density. These unphysical fluctuations 
can be eliminated by mapping the Gaussian random variable onto a positive definite negative 
binomial distribution (NBD). There is a long history of modeling fluctuations in pp collisions 
with negative binomial distributions (for recent examples see EH)). Writing 

m + n) n n k k 

NBD(n, k; n) = T , (2) 

V ; Y(k)Y(n + 1) (n + k) n+k 

where n is the sampled value, h is its mean and k controls its variance, we identify njn with 
e(r)/eo. In the limit of large h (we take n- 100), the NBD becomes a continuous function 
P NBD (;y) of the reduced variable y-njn whose width parameter k we adjust such that its variance 

)= \ + \ agrees with the squared Gaussian width Cov[e(r)/e ]|,=o = (Ae(0)/e ) 2 = (0.684) 2 . 
The mapping is now achieved by replacing each value of x = e/eo from the Gaussian random field 
by a new value y = e new /eo sucn that the cumulative Gaussian distribution function at x coincides 
with the cumulative NBD distribution at njn — e Dew /eo -y: 

f P QwS s(x')dx' = f P mD (y')dy', (3) 

J-co Jo 

The resulting negative binomial random field e(r)/eo is positive definite and retains the two- 
point covariance embedded in the original GRF. A small section of this field and its two-point 
correlation function are shown in Fig. [2] 
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Figure 2: (Color online) Same as Fig.^ but for the NBD random field. Note that in the left panel the plotted quantity 
e(r)/eo _ l is now > — 1 everywhere. 



To imprint these energy density fluctuations onto the initial conditions for a heavy-ion colli- 
sion, we take the transverse energy densities dE/d 2 r ± dy for events generated with the MC-KLN 
model [ 10 1 and multiply them with the NBD random e(r)/eo field taken from randomly selected 
and appropriately sized sections of the final 4000 x 4000 field grid: 

dEfiucXr) = dE KLN (r) ^ e(r) 
d 2 rdy d 2 rdy £q 



3. Results and Conclusions 

This texturing procedure was applied to 20, 000 MC-KLN Au-Au events at V* = 200 A GeV 
(using the MC-KLN code with Gaussian nucleons of width <x = 0.54 fm), partitioned into equally 
sized bins in the number of participants of width AN part = 1000. For both the textured and untex- 
tured events we compute the harmonic eccentricity coefficients e„ using the definition 

[rdrd^r 2 e^ d -^f 

6 " Jrdrd,r 2 ^ 

and average them over the event ensemble. The ratios of these averages are plotted in Fig. [3] for 
the harmonics n = 2, . . . , 5 as functions of A^ part . 

For the nominal correlation length a = 0.28 fm , shown in the left panel, we see that the 
gluon-field fluctuations induce only a small increase in the eccentricity harmonics e„, reaching 
5-10% in central collisions and falling off in more peripheral ones. We should note, however, 
that our fluctuation texture assumes constant (i.e. position-independent) Q s or a, with a value 
expected (on average) in the center of central Au-Au collisions. More realistically, a 2 should 
vary inversely with the position-dependent nuclear thickness function (which, on average, is 
largest in the fireball center): 

a 2 * l/fijoe l/r(r). (6) 

Consequently, the Miiller-Schafer texture with constant a = 0.28 fm provides an approximate 
lower bound on the eccentricity enhancement caused by sub-nucleonic fluctuations in 200 A GeV 
Au-Au collisions. 
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Figure 3: (Color online) Centrality dependence of the gluon-field fluctuation induced enhancement of e„ in 200AGeV 
Au-Au events, for a correlation length ci = 0.28fm as used in |1| (left panel) and a value a = 0.54fm as large as the 
nucleon size (right panel). 



To obtain an estimated upper bound on the eccentricity enhancement caused by sub-nucleonic 
fluctuations, we inflate the correlation length a to the radius of our Gaussian nucleons, cr = 0.54 fm 
and repeat the texturing procedure. (This amounts to reducing Q s by about a factor 2.) As seen 
in the right panel of Fig. [3] with the larger correlation length the gluon fluctuations increase the 
eccentricities e„ by larger factors, reaching now 20-25% in central collisions (and again falling 
off in peripheral ones). For both values of the correlation length, the sub-nucleonic fluctuation 
effects on e„ are strongest in central collisions; in more peripheral collisions, fluctuations in the 
nucleon positions dominate the fluctuation effects on e„. 

To summarize, we generated a toy model to analyze the effects of sub-nucleonic color fluc- 
tuations on the centrality dependent eccentricty harmonics e„ for MC-KLN initial conditions. 
While we qualitatively confirm earlier findings J7J that such fluctuations tend to increase the 
€„, only relatively small enhancements (smaller than those reported in [8]) are found for realistic 
values of the gluon field correlation length. Correlations over larger distances generate larger 
eccentricities. Due to the assumption of a position-independent saturation momentum Q s or cor- 
relation length a, our implementation of sub-nucleonic fluctuations is much less realistic than the 
one in the IP-Glasma model of Refs. EE); it has, however, the advantage of allowing us to turn 
the sub-nucleonic fluctuations on and off at will and thus to study their effects on e„ in isolation. 
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